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' Abstract. The first purpose of tfiis paper is to present a systemic study of some fami- 

, lies of multiple ^-Bernoulli numbers and polynomials by using multivariate q-Volkenborn 

' integral (= p-adic (/-integral) on Zp . From the studies of these (/-Bernoulli numbers and 
polynomials of higher order we derive some interesting (/-analogs of Stirling number 

r identities. 

H 
^' 

. 81. Introduction 



Let q be regarded as either a complex number G C or a p-adic number q G 
, ^ ! Cp. If G C, tiien we always assume IqI < 1. If g G Cp, we normally assume 

, |1 — < p p-i , which implies that q^ = exp(a;logg) for \x\p < 1. Here, | ■ |p is the 

^ i p-adic absolute value in Cp with \p\p = K The g-basic natural number are defined 

(3 . hj [n]q = = 1 + q + ■ ■ ■ + q^~^, ( n G N), and q'-factorial are also defined as 

^ . i^Tijg! = [njg- [n — l]g ■ ■ ■ [2]g- [l]g. In this paper we use the notation of Gaussian binomial 

Q ! coefficient as follows: 

n\ _ [n]q\ _ [n]q ■ [n - l]q- ■ -[n - k + l]q 



X- (1) 



kjq [n-k]ql[k]ql [k]q\ 



Note that lim^^i (^)^ = (^) = ''^'^ ^'i''"' '^'^^K The Gaussian coefficient satisfies the 
following recursion formula: 
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From thus recursion formula we derive 



(3) Q= Yl g'^^+'^^+'-'+^^N see [15, 20, 21]. 

Let p be a fixed prime. Throughout this paper Zp, Qp, C, and Cp will, respectively, 
denote the ring of p-adic rational integers, the field of p-adic rational numbers, the 
complex number field, and the completion of algebraic closure of Qp. For d a fixed 
positive integer (p, d) = 1, let 



X = Xd = limZ/dp^Z, and Xi = Zp, 
X*= U (a + dpZp), 

0<_a<_dp 
(a,p) = l 

a + dp^Zp = {x e X\x = a (mod dp^)}, 



where a e Z lies in < a < dp , cf. [8-18]. For x e Cp, we use the notation 
[x], = ijEj , cf. [1-6]. 

We say that / is a uniformly diflFerentiable function at a point a E Zp and denote 

f(x) - f(y) 

this property by / e UD(Zp), if the difference quotients Ff(x,y) = have 

x-y 

a limit / = /'(a) as {x, y) (a, a). For / G UD{Zp), let us start with the expression 



^ 0<j<p^ 0<i<p^ 

representing a g-analogue of Riemann sums for /, cf. [8, 21-23]. The integral of / on 
Zp will be defined as limit {n oo) of those sums, when it exists. The g-Volkenborn 
integral (=p-adic g-integral) of the function / e UD(Zp) is defined by 

(4) W) = / f{x)dixq{x) = Yua^ E see [8] . 

The Carlitz's g-BernouUi numbers (3h^q can be determined inductively by 

r 1 if = 1 
/3o., = l. <,W+1)'-A., = |„ 

with the usual convention of replacing by Pi^q{see [2, 3, 24, 25]). 

In [8], it was shown that the Carlitz's g-BernouUi numbers can be represented by 
p-adic g-integral on Zp as follows: 

(5) / [x]'^dfXq{x) = / [X]'^dllq{x) = Prn,q, m e Z+. 
JZp Jx 
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The A;-th order factorial of the q'-number [x]q, which is defined by 



is called qf-factorial of x of order k, cf.[15, 19, 20]. Prom this we note that (^)^ = -^^pf? 
cf. [4-16]. The theory of g-number and the factorial of g-number are applicable in 
the many areas related to mathematics , mathematical physics and probability. For 
example, we consider a sequence of Bernoulli trials and assume that the conditionally 
probability of success at the n-th trial, given that k successes occur before that trial 
varies geometrically with n and k. specifically, suppose that the probability of success 
at the n + 1-th trial, given that k success occur up the n-th trial, is given by 

Xn,k = g«"+''^+'^(e M), = 0, 1, 2, • • • , n, n = 0, 1, 2, • • • , 

with a, b and c such that < A^^^ < 1. The particular case 6 = corresponds to the 
assumption that the probability of success at any trial depends only on the number of 
previous trials, while the other particular case a = corresponds to the assumption 
that the probability of success at only trial depends only on the number of previous 
success. The purpose of this paper is to present a systemic study of some families 
of multiple g-BernouUi numbers and polynomials by using multivariate g-Volkenborn 
integral(=]3-adic g-integral) on Zp. From the studies of these g-BernouUi numbers and 
polynomials we derive some interesting qf-analogs of Stirling number identities. That 
is, the g-analogs of many classical Stirling number identities are formulated and their 
interesting features are revealed in this paper. 



§2. Q'-Bernoulli numbers associated with qf-Stirling number identities 

— ^ 

In this section we assume that g e Cp with |1 — q\p < p p-i . From the definition 
of [x]q we can easily derive the following equation. 

qn qx lg*-l 1 

(6) q^'lx - n]q = = [x]g - [n]q, and [-x]g = ^y^^ = "^W^- 

Let {Eh){x) = h{x + 1) be the shift operator. Then we consider the g-difference 
operator as follows: 

n 

(7) = (£; - q'-^l) , where {Ih){x) = h{x). 

i=l 

From (6) and (7), we note that 

(8) /(^) = y. f!! A^/(o), 



1- 

n>0 ^'"^ 1 



where 

(9) ^^m = J2(l) {-^)'l^'^h{n-k). 



The qf-Stirhng number of the second kind is defined by Carhtz as follows: 

■®E(-i)V«(f 

By (9) and (10) we easily see that 



-(2) 

(11) .2(n,/c,g) = ^A,^0^ 
Prom (11) we can also derive the following equation. 

n / \ n (2) —("2'') 

(12) [^r = EU W,!.2(fe, n-K <i)<f^ - Y>\.q ^ C- k\ ! ^ 

fc=0 ^ ' 1 k=0 

By (2), we easily see that 



(13) / ( 1 dfi,{x) = T^fV^ 



^ (-1)" ^(n+l)-("r) 



From (5), (12) and (13) we can derive the following theorem. 
Theorem 1. For m e Z+, we have 



(14) Pm,q = q E Ti:lTr(-l)'^2(A:, n - A;, q), 



k=0 

where Prn,q O'f^ m-th Carlitz q-Bernoulli numbers. 
The qf-Stirling numbers of the first kind is defined as 



(15) (1 - qr[x]n,, = - q--^+'q'-') = E ( J Q^'^h'^YQ 

i=l 1=0 ^ ^ 1 



It is easy to see that 



(16) q'- = ([x],(g - 1) + 1)^ = E - 

m=0 ^ ^ 



im 



From (16) we note that 



By (12), (15) and (17) we obtain the following theorem. 
Theorem 2. For n e Z+ we have 

In [3] Carlitz has given the following relation. 
(18) s,(„,k.,) = to-l)-'=E(-l)'-'(*!")(')")^. 



and 



1 i=o 

By simple calculation we easily see that 



(19) g" =J2(1- l)'?*'' (t) W'.« = E f E (9 - 1)' (I) 9) 

A;=0 ^ / g m=0 \fc=m ^ ^ Q 

By using p-adic qf-integral on Zp we have 

(20) / q-'df.,{t)=j2(l)i<l-'^rPrn,,. 

J^v m=0 

From (19) and (20) we derive 

(21) (l) = E(^-i)-"^^'(^)^^i(^'-'^)- 
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Prom the definition of the first kind Stirhng number we note that 
(22) 5© f^) [n],\ = [a;]„,,9(^) = si(n, /c, J. 

By (13) and (22) we have 

-1 -1 

By (14), (15), (17) and (23) we obtain the following theorem. 
Theorem 3. Forn,j G we have 

'-j 

k=j 

Moreover, 



-1 -1 ^ 



§3. Multivariate p-adic g-integral on Zp associated with g-Stirling numbers 

— ^ 

In this section we also assume that q G Cp with |1 — g|p < p p-i . For any positive 
integers k, m, we consider the following multivariate p-adic qf-integral on Zp related 
to g-BernouUi polynomials of higher order as follows: 
(25) 

(k) / \ (fc) 

In the special — 0, (3n,q{0) — Pn,q wiU be Called the g-BernouUi numbers of 

order k. From (25) we note that 

Thus, we obtain the following theorem. 
Theorem 4. For m,k E Z_|_, we have 

(-1)^0 m k\ 



/^n.U (l_,)nZ. 
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Now we also define Pn,q^ {x) as follows: 



1 

(27) ^'n-:'H^) = TT—^i: 



where n, k are positive integers. Prom (27) we note that 
It is easy to see that 



Q' „ix 



.2q^ (') ^ (fe + n)---(fc + l)fc...(fc-j + l) ^ (t-;) 

^ ' (^■+")n! (j + n)!(/c + n)---(A;+l) C^^^jnl' 

By (27), (28) and (29) we obtain the following theorem. 
Theorem 5. For n, /c e we have 

(30) /'L''(^)=(T^|:(-irCr),f|^." 

From (18) and (30) we derive 



That is, 



1 ^ (-i)^a N,! 



Thus, we note that 



(31) /3(7/)(0) = / • • • / g^ti(fc-^)-.rf/.,(xi) • • -(///.(xfc) 



k\ 



By the same method we see that /^{^^(O) = • • • , Po,q^\Q) = Thus, we 

have 

..(M,,)=[|j/3j-'(0)=||Ll|! = i. 



(k) 

From the definition of /3m,q we can also derive the following equality. 



rr) ki 



Therefore we obtain the following: 
Finally, we observe that 

9^^^ [X]n,q = [xUq[x - l]q • • • q''~'^[x - n + l]g = [x]q ■ {[X]q -!)••• {[X]q - [u - 

Thus, we have 

§4. Further Remarks and Observations 

In this section, let p be a fixed odd prime number. For n e N, /c e Z_(_, we consider 
the following g-Euler numbers of higher order. 

(32) E'^^\x,q)=l ■■■[ [f]x, + x]Jg^"--^(-^)d/x_,(xi)---d/x_,(xJ, 
where H-q{x + p^Zp) = -^^{-qY = [^^' [l^]- (32) we note that 



q)[xt. 



l+q 

k 



[2]" ( — 1)'^'^ 

The g-binomial formulae are known as 

n n y \ 

(33) n(«+^^^"')=E(j 

i=l /c=0 ^ 



and 



n(a-v-T'=E("^r') 

i=l fc=0 ^ ^ 1 

By (32) and (33) we obtain the following: 



Proposition 6. For n e N, A; e Z+, lye have 



Seeking to define a suitable polynomial analogue for negative value of n, we give 
the definition as follows: 

(-) 4-"'(...)^(,^4|:(:)(-i)v^(n(i..-)). 

From (33) and (34) we can also derive the following Eq.(35). 

(35) = ^ J-g (^■)(-l)',-g (;')^,(»,c«).. 
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